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1.0  INTRODUCTION 

Lasers  are  being  used  in  environments  that  may  have  adverse  transmission  properties. 
Communications  applications,  for  example,  may  require  that  the  laser  signal  be  transmitted 
reliably  over  long  distances  that  may  occasionally  be  obscured  by  clouds  of  rain,  dust,  or 
fog. 

One  important  role  of  laser  testing  is  to  define  the  transmission  degradation  that  adverse 
environments  produce.  It  is  essential  that  realistic,  well-defined  test  conditions  be  available 
for  such  tests.  The  purpose  of  the  present  test  effort  was  to  make  transmission 
measurements  for  carbon  dioxide  (CO^)  and  helium-neon  (HeNe)  laser  beams  attenuated  by 
water  drops.  An  existing  AEDC  facility  and  technology  developed  for  other  applications 
were  used  to  obtain  the  measurements.  The  approach  was  to  introduce  attenuating  particles 
into  a laser  beam  and  directly  measure  the  resulting  transmission  change. 

2.0  APPARATUS 


2.1  TEST  UNIT 

The  test  facility  is  a 60-ft  long,  8-ft-diam  steel  tank  that  can  be  evacuated  to  any  desired 
pressure  between  an  atmosphere  and  a few  millimeters  of  mercury  (mmHg).  The  normal  test 
gas  is  air,  but  other  gases  are  available.  Humidity  control  is  also  available,  since  both  a dry 
air  source  and  a steam  line  are  installed.  Liquid-nitrogen  plumbing  is  installed  for  local 
cooling  inside  the  tank  when  required.  The  experimental  arrangement  is  depicted  in  Fig.  1, 
and  the  general  optical  arrangement  is  shown  in  Fig.  2. 

A water  drop  generator  provided  the  particulate  field  for  the  attenuation  experiments. 
Drop  size  was  varied  between  0.3  and  1.4  mm.  A monodisperse  water  drop  generator  was 
mounted  below  the  coincident  laser  beams  so  that  the  trajectory  of  the  stream  of  water  drops 
entered  the  beam  and  then  fell  out  of  the  beam  as  depicted  in  Fig.  3.  The  monodisperse 
generator  used  is  an  example  of  the  equipment  developed  for  hypervelocity  erosion  testing  in 
a ballistics  range.  In  that  application,  the  generator  is  mounted  so  that  the  drops  fall 
vertically  into  the  path  of  the  projectile.  For  laser  attenuation  testing,  the  usual  vertical 
mounting  would  not  place  enough  drops  in  the  beam  to  produce  appreciable  attenuation. 

The  monodisperse  generator  employs  a vibrating  orifice  that  breaks  the  .water  stream 
into  uniformly  sized  drops.  The  orifice  is  vibrated  by  a piezo-electric  crystal,  to  produce  a 
single  drop  for  each  electrical  pulse  from  a signal  generator.  Thus,  the  easily  measured  and 
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Figure  1.  Experimental  layout. 
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controlled  signal  generator  frequency  is  the  drop  formation  frequency.  The  water  drops 
have  a diameter  uniformity  of  ± 2 percent,  and  the  particle  diameter  can  be  calculated  from 
the  operating  conditions  of  the  generator  to  better  than  2-percent  accuracy.  This 
experimental  arrangement  places  a stream  of  water  drops  into  the  coincident  laser  beams,  as 
depicted  in  Fig.  4,  so  that  the  attenuation  of  each  beam  by  a common  collection  of  water 
drops  can  be  measured . 

2.2  INSTRUMENTATION 

A CO2  laser  and  an  HeNe  laser  were  used  for  the  attenuation  experiments.  The  C02  laser 
was  a 3-w  continuous  device.  An  attenuator  was  installed  on  the  laser  to  reduce  the  power 
output  to  the  1/2-w  level.  The  initial  beam  was  expanded  to  a 1.7-in.  diameter  using  a 
germanium  lens  and  was  then  focused  onto  the  detector  at  the  far  end  of  the  lank.  This 
focusing  was  accomplished  by  a 6-in.-diam,  48-in. -focal-length,  spherical,  concave  mirror. 
The  beam  size  at  the  detector  was  approximately  1/4  in.  in  diameter. 

The  HeNe  laser  was  a 15-mw  continuous  device.  The  beam  was  expanded  to  a 1.5-in. 
diameter  and  then  focused  onto  a detector  at  the  far  end  of  the  tank.  The  expanding  and 
focusing  were  accomplished  using  a commerical  beam  expander  attached  directly  to  the 
laser. 

The  geometric  arrangement  is  depicted  in  Fig.  1.  A more  detailed  diagram  of  the  optical 
arrangement  is  given  in  Fig.  2,  in  which  the  angle  between  the  two  laser  beams  is  greatly 
exaggerated.  An  important  feature  of  the  arrangement  is  that  the  two  beams  are  essentially 
coincident  through  the  particulate  region. 

The  output  of  each  detector  was  directed  to  a minicomputer  so  that  the  analog  output 
signals  could  be  digitized  online  and  stored  in  digital  form  for  subsequent  computer 
processing.  The  system  has  the  capability  to  digitize  and  record  both  channels  for  0.32  sec  in 
real  time  at  a frequency  of  3.135  kHz. 

1 

Other  features  of  the  experimental  arrangement  included  a provision  for  completely 
shielding  the  detectors  from  the  laser  beams  to  establish  the  zero  level  output  and  a provision 
for  blocking  the  water  drops  injected  into  the  laser  beams  so  that  the  full-scale  output  of  the 
detectors  could  be  established. 
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Figure  4.  Particulate  region  as  viewed  along  laser  beam  axis. 
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3.0  PROCEDURE 

3.1  EXPERIMENTAL  PROCEDURE 

Test  measurements  consisted  of  direct  measurements  of  the  laser  attenuation  produced 
by  the  water  drops  and  measurements  needed  to  define  the  size  and  population  of 
attenuating  water  drops  in  the  laser  beams.  For  each  test  measurement,  the  water  flow  was 
first  established.  The  water  flow  rate  to  the  drop  generator  was  measured,  as  was  the 
frequency  of  the  signal  supplied  to  the  drop  generator.  These  quantities  gave  the  drop 
formation  frequency  and  size.  Other  measurements  included  the  height  of  the  drop 
trajectory  into  the  laser  beam  (h  in  Fig.  3)  and  length  of  the  drop  trajectory  in  the  laser  beam 
(L  in  Fig.  3).  These  quantities,  along  with  the  drop  formation  frequency  and  the  drop  size, 
defined  the  drop  spacing  along  the  trajectory.  Some  variation  existed  in  the  trajectory  length 
along  the  axis,  L.  This  quantity  was  also  needed  to  define  the  cross-sectional  area  the  drops 
present  to  the  laser  beams. 

After  the  flow  parameters  were  measured,  the  detector  outputs,  which  correspond  to  no 
laser  signal,  were  established  by  completely  shielding  the  detectors  from  the  laser  beams. 
This  information  was  recorded  over  a 0.32-sec  time  span,  so  that  the  effects  of  various 
sources  of  noise,  if  they  existed,  could  be  averaged  out.  Once  the  null-level  output  was 
established,  the  full-scale  output  was  measured  in  a like  manner,  by  blocking  the  water  flow 
into  the  laser  beam.  An  example  of  the  full-scale  output  recorded  for  0.32  sec  is  shown  in 
Fig.  5. 

A test  shot  was  then  made  during  which  the  outputs  of  the  respective  detectors  were 
recorded  for  0.32  sec.  The  test  outputs,  corresponding,  to  the  calibration  values  shown  in 
Fig.  5,  are  shown  in  Fig.  6.  Figures  5 and  6 show  that  the  signal  from  the  HeNe  detector  was 
relatively  uniform,  but  the  signal  from  the  COj  detector  had  a 120-Hz  noise  component, 
which  apparently  was  inherent  in  the  laser.  The  detector  outputs  were  plotted  for  each  test  as 
shown  in  Figs.  5 and  6.  If  they  appeared  to  be  satisfactory,  the  data  points  were  recorded  on 
magnetic  tape  for  subsequent  computer  processing. 

3.2  EXPERIMENTAL  DATA  REDUCTION 

The  transmission  is  defined  as  the  ratio  of  the  power  transmitted  through  the  attenuation 
region  to  that  incident  on  the  attenuation  region.  The  experimental  transmission  is  therefore 
given  by 


T = (V  - V0)  / (VF  - V 0) 


(1) 
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Point  Number 


Figure  5.  Sample  of  full-scale  output  of  detectors. 


Point  Number 

Figure  6.  Sample  output  of  detectors  during  test  shot. 
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where 


V is  the  detector  output  voltage  during  the  test, 

V0  is  the  detector  output  voltage  with  laser  beam  blocked,  and 
Vp  is  the  detector  output  voltage  with  no  attenuation 

The  instantaneous  attenuation  data,  read  from  the  magnetic  tape,  were  averaged  over 
the  recording  interval  of  0.32  sec  to  eliminate  the  effect  of  random  signal  noise.  These 
averaged  values  were  used  in  Eq.  (1)  to  compute  experimental  values  of  transmission. 

3.3  CHARACTERIZATION  OF  ATTENUATING  MEDIUM 


The  characterization  of  the  water  drops  in  the  attenuation  region  has  been  treated  in  an 
approximate  manner.  The  number  in  the  region,  their  size,  and  approximate  geometric 
locations  can  be  determined.  A major  uncertainty  is  the  appropriate  attenuation  area  for 
adjacent  drops  with  overlapping  cross-sectional  areas. 


The  transmission  through  an  attenuation  region  can  be  expressed  as 


T = 


(2) 


where  AL  is  the  cross-sectional  area  of  the  laser  beam  and  AE  is  an  extinction  area 
attributable  to  the  attenuating  effects  of  the  medium  on  the  laser  beam.  The  extinction  area 
is  related  to  the  cross-sectional  area  the  water  drops  present  to  the  beam  by  an  extinction 
efficiency,  QE, 


At  - 0E  E A, 


1-1 


(3) 


where  A,  is  the  effective  cross-sectional  area  of  the  ith  particle  and  N is  the  total  number  of 
particles  in  the  laser  beam.  The  problem  is  to  determine  the  individual  drop  contributions, 
A,,  from  the  trajectory  measurements,  flow  rate,  and  drop  formation  frequency. 


The  time  interval  between  drops  is  given  by 


where  f is  the  drop  formation  frequency.  The  drop  diameter,  d,  is  given  by 


(4) 


(5) 


\ 
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with  the  volumetric  flow  rate  denoted  by  Q.  The  number  of  drops  in  the  attenuation  region 
is  given  by 


N = 2f  JljT  (6) 

v B 

with  the  maximum  drop  trajectory  height  into  the  beam  given  by  h. 

For  the  special  case  in  which  the  drops  are  well  separated  so  that  no  interactions  need  be 
considered,  the  effective  cross-sectional  area  would  be  given  by 


In  the  present  case,  overlapping  area  between  adjacent  drops  must  be  taken  into  account. 
When  all  drops  follow  a common  trajectory,  the  vertical  spacing  between  adjacent  drop 
centers  is  given  by 


'd  [Vv  - 3 i2t  - td> 


(8) 


with  the  initial  vertical  velocity  given  by  Vv  and  the  time  after  injection  of  the  drop  into  the 
beam  given  by  t.  For  a common  trajectory,  the  extent  of  the  overlapping  areas  between 
adjacent  drops  can  be  determined  by  Eq.  (8).  The  result  of  the  calculation  of  the  overlap 
area  between  adjacent  drops  is  shown  in  Fig.  7. 


The  horizontal  extent  of  the  drop  trajectory  in  the  laser  beam  is  denoted  by  L in  Fig.  3. 
Experimentally,  for  each  test  a certain  amount  of  jitter  was  found;  therefore,  the  dimension 
L is  an  average  distance  over  which  the  drops  remain  in  the  laser  beam.  The  jitter  in  this 
dimension  was  used  to  estimate  the  magnitude  of  the  random  drop-to-drop  dispersion  from 
the  mean  trajectory.  This  scatter  must  be  considered  in  the  analysis  procedure,  since  it  enters 
into  the  calculation  of  the  overlapping  area  between  adjacent  drops. 


The  jitter  is  assumed  to  be  caused  a random  variation  in  the  initial  velocity  of  the  water 
drops  entering  the  laser  beams.  Under  this  assumption,  the  fractional  variation  in  L is  twice 
the  fractional  variation  in  the  initial  velocity,  and  the  random  component  in  the  vertical  drop 
position  is  expressed  as 


by  - 


^ G5L.) 

t 


(9) 
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Figure  7.  Calculation  of  overlapping  projected  area  of  adjacent  drops. 

where  SL  is  the  extent  of  the  variation  in  L,  Vv  is  the  mean  initial  vertical  drop  velocity, 
and  t is  the  time-of-flight  of  the  drop. 

A comparable  spread  in  the  horizontal  location  of  the  drops  perpendicular  to  the  laser 
beams  is  also  present.  Equation  (9)  also  estimates  the  extent  of  the  horizontal  jitter  of  the 
individual  drops  with  respect  to  the  vertical  plane  of  the  main  trajectory. 

The  projected  cross-sectional  area  of  the  water  drops  has  been  calculated  for  each  test  in 
the  following  manner.  The  first  drop  is  assumed  to  be  centered  on  the  border  of  the  laser 
beam.  Half  of  its  area  is  counted  as  being  in  the  laser  beam.  The  vertical  location  of  the  next 
drop  is  then  calculated  by  assuming  that  it  moves  along  the  mean  trajectory  for  a drop 
formation  time.  Equation  (9)  then  establishes  a maximum  deviation  from  the  mean 
trajectory  likely  to  occur.  A uniformly  distributed  random  number  between  0 and  5y  as 
given  by  Eq.  (9)  is  taken  to  be  the  magnitude  of  the  deviation  in  position  from  the  mean 
trajectory.  Another  uniformly  distributed  random  number  between  0 and  2ir  denotes  the 
direction  of  the  deviation  from  the  mean  trajectory.  This  quantity  is  added  vectorally  to  the 
vertical  spacing  calculated  for  the  mean  trajectory  to  give  the  separation  of  the  adjacent 
drops.  The  area  of  the  second  drop  is  added  to  the  contribution  of  the  first  and  the  overlap 
area  (if  overlap  exists)  is  subtracted  from  the  accumulated  area.  The  overlap  area  is 
calculated  as  shown  in  Fig.  7,  with  the  calculated  drop  separation  distance  used  as  hd. 
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The  contribution  of  the  next  drop  to  the  accumulated  area  is  then  considered  in  a like 
manner,  until  the  last  drop  considered  falls  out  of  the  laser  beam.  This  calculation  gives  the 
total  number  of  drops  in  the  laser  beams  and  the  appropriate  total  cross-sectional  area 
occupied  by  the  water  drops. 

3.4  TRANSMISSION  PREDICTIONS 

The  experimental  results  have  been  analyzed  using  the  simplified  approximate  approach 
presented  in  Ref.  1.  A sphere  placed  in  a light  beam  is  assumed  to  interact  with  the  incident 
radiation  in  the  same  manner  as  an  opaque  disk  of  the  same  diameter.  The  first  effect  is  that . 
all  the  radiation  intercepted  by  the  disk  is  blocked.  By  Babinet’s  principle  (Ref.  2)  a like 
amount  of  radiation  is  diffracted  from  the  rim  of  the  disk.  Thus,  the  incident  radiation  on  a 
circular  region  of  diameter  V2d  is  scattered  from  its  original  path  by  the  presence  of  a disk  or 
sphere  of  diameter  d. 

For  absorbing  spheres  with  sufficiently  large  size  parameters,  i.e.  (Ref.  3), 


x = rd  A > 2000 

the  scattered  radiation  intensity  distribution  can  be  analyzed  using  diffraction  theory. 
However,  for  small  scattering  angles,  diffraction  theory  appears  to  approach  the  results  of 
the  much  more  elaborate  Mie  scattering  theory  for  values  of  x as  small  as  80  (Ref.  4).  For  the 
present  experiments,  80  < x < 4000,  with  scattering  angles  of  interest  of  less  than  a degree. 
Therefore,  diffraction  theory  should  be  appropriate  to  the  present  experimental  conditions. 

The  power  received  by  the  detector  with  no  intervening  water  drops  is 


P 


LA, 


(10) 


with  the  incident  radiation  intensity  at  the  attenualion  region  given  by  I0  and  the  cross- 
sectional  area  of  the  laser  beams  at  the  region  denoted  by  Al.  When  a water  drop  occupies 
the  attenuation  region,  the  power  received  by  the  detector  is 


P 


'o'  I-  ” W 0A ; 


<?• 


(ID 


where  Ap  is  the  cross-sectional  area  of  a water  drop,  Ad  is  the  cross-sectional  area  of  the 
detector,  and  ls  is  the  scattered  radiation  intensity.  The  scattering  efficiency  of  the  drop  is 
denoted  by  Q^.  The  first  contribution  to  the  power  received  by  the  detector  in  Eq.  (11)  is  that 
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attributable  to  the  power  incident  on  the  attenuation  region.  The  second  term  represents  the 
power  scattered  by  the  water  drop,  and  the  third  term  represents  the  scattered  power  that  is 
diverted  toward  the  detector, 

The  scattering  efficiency,  Qs,  has  been  evaluated  using  Eq.  (12),  which  has  been  taken 
from  Ref.  4: 


Qs 


= 2 - 


4 4 

— sinpH (1  — uos  p) 

P P 2 


02) 


where 


P = 1)  (13) 

The  particle  diameter  is  d,  the  radiation  wavelength  is  X,  and  m is  the  refractive  index  of  the 
water  drop  (taken  to  be  1.33). 

Equations  (12)  and  (13)  have  been  evaluated  for  the  range  of  the  test  parameters.  The 
results  are  shown  in  Fig.  8 as  the  envelope  of  a rapidly  oscillating  function  of  particle 
diameter. 

The  transmission,  T,  is  given  by  Eq.  (14): 


The  power  scattered  in  the  direction  of  the  detector  can  be  evaluated  using  Fraunhofer 
diffraction  theory  (Ref.  3): 


(15) 


where  i(x,0)  is  the  Fraunhofer  intensity  function,  0 is  the  scattering  angle,  f is  the  distance 
from  the  scattering  particle  to  the  detector,  and 


x = n d/  it 


(16) 
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Figure  8.  Envelope  of  scattering  efficiencies  calculated  from  Eq.  (12). 


The  intensity  function,  i(x,0)  is  given  by  (Ref,  2) 


ifv.  9) 


(x  Sill  9) 


"|2 


in  9 


(17) 


Values  of  the  intensity  function  for  the  region  of  interest  are  tabulated  in  Tables  1 and  2. 
The  zeros  of  the  intensity  function  occur  at  values  of  x sin  6,  which  correspond  to  the  roots 
of  J]»  the  first  order  Bessel  function.  For  experiments  at  either  laser  wavelength,  the  angular 
extent  of  the  detector  is  less  than  the  scattering  angle  corresponding  to  the  first  root  of  Jj. 

The  intensity  function  must  be  integrated  over  the  area  subtended  by  the  detector  to 
evaluate  the  fraction  of  the  scattered  radiation  that  is  scattered  in  the  direction  of  the 
detector.  Using  ( for  the  distance  from  the  attenuation  region  to  the  detector,  dd  for  the 
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detector  diameter,  and  dL  for  the  laser  beam  diameter,  the  integral  of  the  intensity  function 
can  be  expressed  as 


The  results  from  Eq.  (18)  are  shown  in  Fig.  9 for  both  wavelengths  and  the  range  of 
particle  diameters  covered  by  the  present  tests.  The  results  show  that  for  a wavelength  of 
10.6  pm,  very  little  radiation  is  scattered  toward  the  detector,  whereas  for  the  HeNe  laser 
beam  a considerably  larger  fraction  of  the  energy  is  directed  toward  the  detector. 


Table  1.  Fraunhofer  Intensity  Functions,  A = 0.6328  p 


Angle,  deg 


Diameter, 

mm 

_C 

)_ 

0.0065 

0.0130 

0.0195 

0.0260 

0.3 

1 .2302 

X 

1012 

1 ,2214 

X 

io12 

1.1954 

X 

1012 

1.1532 

X 

io12 

1 .0962 

X 

io12 

0.4 

3.8879 

X 

io12 

3.8388 

X 

io12 

3.6947 

X 

io12 

3.4645 

X 

,o12 

3.1623 

X 

io12 

0.5 

9.4919 

X 

io12 

9.3053 

X 

io12 

8.7636 

X 

io12 

7.9191 

X 

io12 

6.8508 

X 

io12 

0.6 

1 .9682 

X 

io13 

1.9127 

X 

io13 

1 .7539 

X 

1013 

1 .5135 

X 

io13 

1.2233 

X 

io13 

0.7 

3 . 6464 

X 

io13 

3.5070 

X 

io13 

3.1150 

X 

io13 

2.5426 

X 

io13 

1.8889 

X 

io13 

0.8 

6.2206 

X 

io13 

5.9115 

X 

io13 

5.0597 

X 

io13 

3.8661 

X 

io13 

2.5908 

X 

io13 

0.9 

9.9642 

X 

io13 

9.3411 

X 

to13 

7.6623 

X 

io13 

5.4194 

X 

io13 

3.2031 

X 

io13 

1.0 

1.5187 

X 

io14 

1.4022 

X 

io14 

1.0961 

X 

io14 

7.0877 

X 

io13 

3.5933 

X 

to13 

1.1 

2.2235 

X 

io14 

2.0185 

X 

io14 

1.4950 

X 

to14 

8.7153 

X 

io13 

3.6558 

X 

io13 

1.2 

3.1492 

X 

1014 

2.8062 

X 

io14 

1 .9572 

X 

io14 

1 .0124 

X 

io14 

3.3452 

X 

io13 

1.3 

4.3376 

X 

io14 

3.7877 

X 

io14 

2.4716 

X 

io14 

1.1135 

X 

io14 

2.6976 

X 

io13 

1.4 

5.8343 

X 

io'4 

4.9841 

X 

io14 

3.0222 

X 

to14 

1.1599 

X 

io14 

1.8354 

X 

io13 

1.5 

7.6884 

X 

io14 

6.4144 

X 

io14 

3.5882 

X 

io14 

1.1420 

X 

io14 

9 . 5035 

X 

io12 
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Diameter, 

imn 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.1 

1.2 

1.3 

1.4 

1.5 


Table  2.  Fraunhofer  Intensity  Functions,  X = 10.6  ju 

Angle,  deg 


-9-  0-0065  0.0130  0.0195  0.0260 


1.5624 

X 

107 

1.5624 

X 

io7 

1.5623 

X 

io7 

1 .5621 

X 

io7 

1.5618 

X 

io7 

4.93B1 

X 

107 

4.9378 

X 

io7 

4.9372 

X 

to7 

4.9360 

X 

io7 

4.9345 

X 

,o7 

1.2056 

X 

108 

1.2055 

X 

io8 

1.2052 

X 

IO8 

1.2048 

X 

io8 

1.2042 

X 

!08 

2.4999 

X 

ID8 

2 . 4996 

X 

io8 

2.4989 

X 

io8 

2.4976 

X 

io8 

2.4958 

X 

io8 

4.6314 

X 

108 

4.6307 

X 

io8 

4.6288 

X 

io8 

4.6256 

X 

108 

4.6211 

X 

io8 

7.9001 

X 

108 

7.8995 

X 

io8 

7.8952 

X 

io8 

7.8880 

X 

10s 

7.8781 

X 

io8 

1 .2656 

X 

,o* 

1.2653 

X 

109 

1.2644 

X 

109 

1.2630 

X 

io9 

1 .2609 

X 

109 

1 .9289 

X 

io9 

1.9284 

X 

109 

1.9267 

X 

io9 

1 .9240 

X 

io9 

1.9202 

X 

io9 

2.B241 

X 

109 

2.8232 

X 

io9 

2.8203 

X 

io9 

2.8155 

X 

io9 

2.8087 

X 

109 

3.9999 

X 

109 

3.9982 

X 

109 

3.9933 

X 

io9 

3.9852 

X 

io9 

3.9739 

X 

109 

5.5092 

X 

109 

5.5066 

X 

109 

5.4987 

X 

109 

5.4856 

X 

io9 

5.4672 

X 

IO9 

7.4102 

X 

io9 

7.4061 

X 

109 

7.3938 

X 

io9 

7.3733 

X 

io9 

7.3447 

X 

io9 

9.7652 

X 

io9 

9.7590 

X 

109 

9.7404 

X 

io9 

9.7094 

X 

109 

9.6663 

X 

IO9 

Figure  9.  Contribution  of  forward-scattered  radiation  to  detector  signal. 
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3.5  UNCERTAINTY  ESTIMATES 


The  independently  measured  quantities  along  with  their  estimated  uncertainties  are  listed 
in  Table  3.  The  assigned  uncertainties  are  based  on  calibration  data  and  experience  and 
represent  an  estimated  standard  deviation  of  the  measured  value  from  the  true  value  of  the 
quantity.  The  independent  quantities  can  be  separated  into  those  that  affect  the  predicted 
transmissions  and  those  that  influence  the  measured  transmissions.  This  separation  is  made 
in  Table  3. 

Table  3.  Estimates  of  Uncertainty  in  Measured  Quantities 


A.  Quantities  Affecting  Predicted  Transmission 

Quantity  Estimated 

Symbol  Description Uncertainty 


h 

Height  of  water  drop  trajectory  in  laser  beams 

0.05 

in. 

L 

Length  of  water  drop  trajectory  in  laser  beams 

0.1 

in. 

Q 

Water  flow  rate 

0.01 

x Q 

f 

Frequency  of  drop  formation 

0.01 

x f 

«L 

Drop-to-drop  variation  from  mean  value  of  L 

0.05 

in . 

B. 

Quantities  Affecting  Measured  Values  of  Transmission 

V 

0 

Voltage  output  of  detector  with  no  incident 
radiation 

0.01 

V 

VF 

Voltage  output  of  detector  with  no  attenuation 
of  Incident  radiation 

0.01 

V 

V 

Voltage  output  of  detector  with  attenuation 
of  incident  radiation 

0.01 

V 

3.5.1  Estimates  of  Uncertainly  in  Predicted  Transmission 


The  magnitude  of  the  uncertainty  in  the  predicted  transmission  can  be  expressed  as 


AT  = 


Ah]5  + [(£)*'!  * [(s)i<3I  * [(I?)"! 


Calculation  of  the  various  terms  for  a number  of  sample  runs  show  that  the  variation  of 
transmission  with  drop  trajectory  length,  L,  is  zero.  The  magnitudes  of  the  other  terms  are 
given  in  Table  4.  For  the  limited  number  of  sample  runs  considered,  the  calculated 
uncertainties  in  predicted  transmission  vary  from  ±0.0309  to  ±0.0723. 
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Table  4.  Estimates  of  Uncertainty  in  Predicted  Transmission  for  Sample  Runs 


Run  No. 

11  Ah 

3h 

11  AQ 

19 

3T  Af 

af 

ii  A4L 

SAL 

AT 

TO 

(-0.26) (0.05) 

(-0.000824X0.607) 

(-0.000137X7.31) 

(-0. 560X0.05) 

±0.0309 

20 

(-0.190) (0.05) 

(0X0.273) 

(0.0000684X7.2) 

(-0.630) (0.05) 

±0.0329 

30 

(-0.120) (0.05) 

(-0.0231) (0.065) 

(0X52.2) 

(-1.44) (0.05) 

10.0723 

40 

(-0.230) (0.05) 

(-0.224) (0.0647) 

(-0.000269X50.2) 

(-0.83) (.05) 

±0.0474 

3.5.2  Estimates  of  Uncertainty  in  Measured  Transmission 


The  measured  transmission  has  been  given  by  Eq.  (1).  The  uncertainty  in  measured 
transmission  is  given  by 


AT 


I r fsT\  wi 2 iy*t  Yw  i 2 r fdj  V\v  i 

yl  L(^/  j + l_(^v  V°J  1 L\^AVfJ 


And  from  Eq.  (1), 


<5T  l 

dv  " v.  -v 

^ 0 


(20) 


(21) 


r}~\  V — \ p 

^ " CVF-tJ1 


(22) 


ii  v.-y  , 

«'V ' ,23) 

Uncertainties  in  measured  transmission  have  been  assessed  using  an  uncertainty  of  0.01 
for  all  voltages  and  Eqs.  (21)  through  (23)  in  Eq.  (20).  The  results  are  tabulated  in  Table  5 
for  the  sample  cases  considered.  For  the  HeNe  laser,  0.00167  < AT  < 0.00261,  and  for  the 
CO,  laser,  the  uncertainty  is  0.00170  < AT  < 0.00469.  These  values  are  small  compared  to 
the  uncertainties  in  predicted  transmission.  The  major  source  of  the  data  uncertainty,  from 
the  uncertainty  analysis,  is  the  characterization  of  the  attenuating  medium.  The  effective 
cross-sectional  area  of  the  water  drops  is  the  least  well-defined  parameter  and  the  drop-to- 
drop  jitter  in  the  trajectory  is  the  major  source  of  this  possible  error. 
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Table  5.  Estimates  of  Uncertainty  in  Measured  Transmission  for  Sample  Runs 

Rur*  Ko . 


10 20 30 40 


IteNe 

co^ 

HeNe 

co2 

HeNe 

C02 

HeNe 

OJl 

8 1 

5T/3V 

0.1317 

0.1572 

0.2014 

0.1678 

0.1661 

0.3421 

0.1696 

0.3106 

3T/3V 

o 

-0.034B9 

-0.02275 

-0.03964 

-0.01934 

-0.02129 

-0.02270 

-0.02631 

-0.02785 

ST/JVp 

-0.09676 

-0.1345 

-0.1618 

-0.01484 

-0.1448 

-0.3194 

-0.1433 

-0.2828 

it 

±0.00167 

±0.00208 

±0.00261 

+0.00170 

±0.00221 

±0.00469 

±0.00224 

±0.00421 

4.0  RESULTS 

The  transmission,  given  by  Eq.  (14),  has  been  calculated  for  each  test  run,  using  Eqs. 
(17)  and  (18)  to  evaluate  the  contribution  of  the  scattering  by  the  particles.  The  results  are 
tabulated  in  Table  6,  along  with  the  experimental  transmission  measurements.  The  drop 
diameter  has  been  calculated  for  each  test  using  Eq.  (5),  and  the  total  number  of  drops  in  the 
laser  beam  has  been  computed  using  Eq.  (6).  The  appropriate  total  projected  area  of  the 
drops  in  the  beam  has  been  calculated  for  each  test  run  assuming  interactions  only  between 
adjacent  drops  as  previously  described. 

For  the  HeNe  laser  Fig.  10  compares  the  predicted  values  obtained  from  Eq.  (14)  with 
the  experimental  transmission  measurements  on  each  test.  Figure  1 1 gives  the  corresponding 
comparisons  for  the  C02  laser.  The  standard  deviation  of  the  measured  data  from  the 
predicted  values  is  also  included  in  Figs.  10  and  11.  Figure  12  compares  the  experimental 
transmission  of  the  HeNe  laser  beam  with  the  experimental  transmission  of  the  C02  laser 
beam  while  both  are  attenuated  by  a common  water  drop  Held.  With  the  exception  of  a 
couple  of  points,  the  agreement  is  within  about  10  percent. 
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Table  6.  Tabulated  Test  Results 


Drop 

Diameter, 

mm 


1 .420 


1.428 


1.426 


1.427 


1.429 


1.435 


1.367 


1.373 


1.373 


1 .383 


1.355 


1,369 


1.023 


1.011 


1.049 


1.053 


1.054 


1.054 


1.038 


1.065 


1.029 


1.041 


1.045 


0.338 


0.326 


Total 
No.  of 
Drops 


Total 
Dr  op 
Area, 
cm2 


99 


0.477 


0.523 


0.452 


0.458 


0.606 


0.587 


0.510 


0.510 


0.516 


0.581 


0.826  j 
0.594  1 


Calculated 

Transmission 


.609 


0.812 


0.794 


0.822 


0.820 


0.761 


0.769 


0.799 


0.799 


0.797 


0.771 


0.611 


0.812 


0.795 


0.822 


0.819 


0.759 


0.768 


0.797 


0.797 


0.797 


0.772 


0.675  0.672 
0.766  0.764 


0.794 


0.490 


0.807 


0.805 


0.805 


r — — ■ ■ — 

Measured 

Transmission 

HeNe 

o 

o 

to 

0.646 

0.601 

0.743 

0.834 

0.781 

0.767 

0.815 

0.806 

0.745 

0.870 

0.782 

0.741 

0.800 

0.794 

0.826 

0.908 

0.818 

0.778 

0.856 

0.751 

0.816 

0.701 

0.772 

0.705 

0,738 

0.735 

0.794 

0.815 

0.866 

0.964 

0.909 

0.746 

0.911 

0.855 

0.733 

0.987 

0.718 

0.807 

0.884 

0.806 

0.343 

0.824 

0.845 

0.732  , 

0.946 

0.800 

0.845 


'0.759 
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Drop 

Diameter, 

mm 


0.328 


0.325 


0.337 


0.340 


0.341 


0.331 


0.350 


0.346 


0.340 


0.334 


0.323 


0.323 


0.326 


0.342 


0.345 


0.350 


0.343 


0.320 


0.336 


0.343 


Total 
No.  of 
Drops 


0.426 


0.142 


0.413 


0.265 


0.1742 


0.361 


0.406 


0.265 


0.129 


0.381 


0.252 


0.174 


0.116 


0.161 


0.290 


0.181 


0.129 


0.310 


0.355 


0.097 


Calculated 

Transmission 

Measured 

Transmission 

HeNe 

co2 

HeNe 

co2 

0.832 

0.837 

0.889 

0.768 

0.944 

0.945 

0.998 

0.976 

0.837 

0.837 

0.778 

0.923 

0.896 

0.894 

0.869 

0,913 

0.931 

0.930 

0.928 

0.879 

0.857 

0.861 

0.907 

0.824 

0.840 

0.837 

0.929 

0.702 

0.896 

0.892 

0.837 

0.951 

0.949 

0.948 

0.949 

0.979 

0.850 

0.852 

0.878 

0.830 

0.901 

0.902 

0.897 

0.884 

0.931 

0.932 

0.901 

0.949 

0.954 

0.955 

0.996 

0.902 

0.936 

0.935 

0.882 

0.896 

0.885 

0.882 

0.904 

0.853 

0.929 

0.927 

0.934 

0.894 

0.949 

0.948 

0.920 

0.985 

0.878 

0.877 

0.863 

0.875 

0.860 

0.861 

0.864 

0.900 

0.962 

0.961 

0.972 

0.936 
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Figure  10.  Comparison  of  measured  and  predicted  transmission,  X = 0.6328  p. 


Figure  11.  Comparison  of  measured  and  predicted  transmission,  X = 10.6  p. 
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Figure  12.  Comparisons  of  transmission  measurements  for  HeNe  and  C02  lasers. 

5.0  SUMMARY  AN1)  CONCLUSIONS 

The  transmission  of  a HeNe  laser  beam  and  a coincident  CO2  laser  beam  through  a 
stream  of  uniformly  sized  water  drops  were  simultaneously  measured.  The  detector  outputs 
were  digitized  online  with  a thousand  data  points  on  each  channel  recorded  for  a duration  of 
0.3  sec.  These  signals  were  averaged  over  the  time  interval  to  remove  any  transient  effects. 
The  online  transmission  data  were  processed  by  computer  to  give  measured  transmission 
data. 

The  attenuation  environment  was  characterized  by  various  measurements  from  which 
the  drop  diameter,  effective  drop  cross-sectional  area,  and  number  of  drops  in  the  laser 
beams  could  be  determined.  Drop  sizes  ranged  from  0.3  mm  up  to  1.4  mm  in  diameter. 

Fraunhofer  diffraction  theory  was  used  to  predict  the  experimental  results.  The 
attenuation  effect  of  a single  water  drop  was  calculated  from  simple  diffraction  theory.  This 
effect  was  converted  into  an  extinction  efficiency,  and  the  extinction  efficiency  was  in  turn 
multipilieti  by  the  total  effective  cross-sectional  area  occupied  by  the  drops  to  yield  a 
predicted  attenuation  for  each  test  run. 
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Comparisons  between  the  measured  and  predicted  transmissions  show  satisfactory 
agreement.  The  magnitude  of  the  deviation  of  the  predicted  transmission  from  the 
experimental  transmission  is  the  same  as  would  be  expected  from  the  uncertainty  in  the 
quantification  of  the  attenuating  medium.  Within  the  limits  imposed  by  the  experimental 
technique,  simple  diffraction  theory  adequately  predicts  the  experimental  results. 

It  is  concluded  that  the  major  source  of  uncertainty  in  laser  attenuation  measurements  is 
the  definition  of  the  attenuating  medium.  To  make  more  definitive  attenuation 
measurements,  a better  defined  attenuating  medium  is  required.  The  ability  to  place  more 
particles  in  the  attenuation  region  would  also  be  desirable. 
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NOMENCLATURE 

Aj  Cross-sectional  area  of  detector 

Ag  Extinction  area 

A,  Effective  cross-sectional  area  of  the  ilh  drop 

Al  Cross-sectional  area  of  laser  beam  at  attenuation  region 

A0  Overlap  area  between  adjacent  drops 

Ap  Cross-sectional  area  of  scattering  panicle 

d Drop  diameter 

(1^  Diameter  of  detector 

dL  Diameter  of  laser  beam  in  attenuation  region 

f Frequency  of  drop  formation 

g Acceleration  of  gravity 

h Height  of  drop  trajectory  in  laser  beam 

hd  Vertical  spacing  between  adjacent  drops 

I0  Power  per  unit  area  incident  on  attenuation  region 

Is  Power  per  unit  area  scattered  by  particle 

i(x,0)  Fraunhofer  intensity  function 
L Length  of  drop  trajectory  in  laser  beam 

m Index  of  refraction  of  particle 

N Total  number  of  water  drops  in  laser  beam 
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P Power  incident  on  detector  ( 

P0  Power  incident  on  detector  with  no  attenuation 

Q Volumetric  flow  rate  of  water  to  drop  generator 

Qe  Extinction  efficiency 

Qs  Scattering  efficiency 

T Transmission 

t Time  counted  from  entry  of  drop  into  laser  beam 

td  Time  internal  between  formation  of  successive  drops 

V Detector  output  voltage  during  test 

VF  Detector  output  voltage  with  no  attenuation 
Vv  Initial  vertical  velocity  of  water  drop 

Vo  Detector  output  voltage  with  laser  beam  blocked 

x Size  parameter  (x  = ird/\) 

5L  Magnitude  of  drop-to-drop  variation  in  trajectory  length  in  laser  beam 

Sy  Random  component  of  vertical  drop  position 

0 Angle  between  laser  beam  axis  and  point  on  detector,  viewed  from  attenuation 

region 

X Laser  wavelength 

q Function  of  size  parameter  and  index  of  refraction  of  particle 
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